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In this paper, we give a parameterization of the class of bivariate symmetric 
orthonormal scaling functions with filter size 6x6 using the standard dilation 
matrix 21. In addition, we give two families of refinable functions which are not 
orthonormal but have associated tight frames. Finally, we show that the class of 
bivariate symmetric scaling functions with filter size 8x8 can not have two or more 
vanishing moments. 
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1. INTRODUCTION 

The most common wavelets used for image processing are the tensor-product of univari- 
ate compactly supported orthonormal wavelets. Of this class of wavelets, only the Haar 
wavelet is symmetric which gives its associated filter the property of linear phase. Since 
Daubechies' work[5], numerous generalizations of wavelets have been developed including 
biorthogonal wavelets, multiwavelets, and bivariate wavelets. Since 1992, several exam- 
ples of bivariate compactly supported orthonormal and biorthogonal wavelets have been 
constructed. See Cohen and Daubechies'93 [3] for nonseparable bidimensional wavelets, 
J. Kovacevic and M. Vetterli'92[8] for nonseparable filters and wavelets based on a gener- 
alized dilation matrix, He and Lai' 97 [6] for the complete solution of bivariate compactly 
supported wavelets with filter size up to 4 x 4, Belogay and Wang'99[2] for a special 
construction of bivariate nonseparable wavelets for any given regularity, and Ayache'99 [1] 
for nonseparable dyadic compactly supported wavelets with arbitrary regularity. See also 
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Cohen and Schlenker'93[4], Riemenschneider and Shen'97[l 1], and He and Lai'98[7] for 
bivariate biorthogonal box spline wavelets. 

It is well-known that in the univariate setting, there does not exist symmetric compactly 
supported orthonormal wavelets except Haar for dilation factor 2. We are interested in the 
construction of symmetric wavelets in the bivariate setting with dilation matrix 21 which 
have compact support and vanishing moments. We start with a scaling funtion </>. Let 

oo 

^ 1 ,c 2 )=n™(e wi/2fc ,^ /2fc ) 

k=l 

be the Fourier transform of </>, where 

N N 
j=0 k=0 

is a trigonometric polynomial satisfying m(l, 1) = 1. In addition, the trigonometric 
polynomial m(x, y) satisfies the orthonormality condition 

\m(x,y)\ 2 + \m(-x,y)\ 2 + \m(x, -y)\ 2 + \m(-x, -y)\ 2 = 1. 

Let ipi(x, y) be the corresponding wavelet satisfying 

Mu u u 2 ) = m t {e^l 2 ,e^l 2 )4>{u 1 /2^ 2 /2),i = 1,2,3, 

where the m, are trigonometric polynomials such that the following matrix 

" m(x,y) m(-x,y) m(x,-y) m(-x,-y) ' 

mi(x,y) mi(—x,y) mi(x,—y) mi(—x,—y) 

m 2 (x,y) m 2 {-x,y) m 2 (x,-y) m 2 (-x,-y) 

_m 3 (x,y) m 3 (-x,y) m 3 (x,-y) m 3 (-x, -y) _ 

is unitary. Moreover, we are interested in symmetric scaling functions cj) with a certain 
number of vanishing moments in the sense that their associated trigonometric polynomial 
m(x, y) satisfies 

m(l/x, l/y) = x~ N y- N m(x,y) 

as well as 

= 0, 0<fc<M-l. 

!/=-! 

The symmetry condition provides (f> with the property of linear phase and the vanishing 
moment conditions provide <f> with polynomial reproduction up to degree M—l. lfm(x,y) 
satisfies the symmetric property, then the associated wavelets can easily be found by using 
<f> and 



dx k 



m(x,y) 



x=-1 



mi(x,y) 
m 2 (x,y) 
m 3 (x,y) 



— m(—x,y) 

— x ■ m(x, —y) 

= x ■ m(—x, —y). 
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In summary, we are looking for trigonometric polynomials m(x, y) = CjkX^y k 

j=0 k=0 

which satisfy the following properties: 

(i) Existence: m(l, 1) = 1 . 

(ii) Orthogonality: \m(x,y)\ 2 + \m(-x, y)\ 2 + \m(x, -y)\ 2 + \m(-x, -y)\ 2 = 1 . 

(iii) Symmetry: m(l/x,l/y) = x~ N y~ N m(x, y). 

(iv) M vanishing moments: m(x, y) = (x + 1) M (y + l) M m(x, y) where m(x, y) is 
another trigonometric polynomial. 

In this paper, we construct a complete parameterization of all trigonometric polynomials 
m(x,y) which satisfy the symmetry condition, the vanishing moment condition, and the 
orthonormality condition for TV = 5 and M = 1. Within this class, we identify a two- 
parameter family which contains the trigonometric polynomials associated with scaling 
functions. Outside of this two-parameter family, we show that the remaining trigonometric 
functions are not associated with scaling functions but instead determine families of tight 
frames. Finally, we show that there are no trigonometric polynomials for TV = 7 and 
M = 2, and consequently no symmetric bivariate scaling functions with two vanishing 
moments for the support size we are considering. 

The paper is organized as follows. Section 2 gives the parameterized solution when 
TV = 5 and M = 1. The problem is broken down into four cases which are dealt with in 
turn. Section 3 discusses the orthonormality of the solutions from Section 2 and concludes 
with a numerical experiment comparing Haar, D4, and one solution from Section 2. The 
last two sections show that these trigonometric polynomials cannot have higher vanishing 
moments(i.e M > 2) even for N = 7. 



2. THE 6x6 CASE 

Our goal is to parameterize the coefficients of the trigonometric polynomials which 
satisfy properties (i)-(iv). We begin our investigation with trigonometric polynomials 
whose filter size is 6 x 6 with one vanishing moment, i.e. N = 5 and M = 1. The case 
N = 1 is trivially the tensor product Haar function, and the case N — 3 has eight singleton 
solutions given by He and Lai'97 [6]. 

Let us express m(x,y) in its polyphase form, i.e., 

m{x,y) = f a (x 2 ,y 2 )+xfb(x 2 ,y 2 )+yf c (x 2 ,y 2 ) + xyf d (x 2 ,y 2 ) 



a 


bo 


ai 




CL2 


b 2 - 




1 " 


CO 


d 


Cl 




C2 


d 2 




X 


«3 


h 


a4 


b A 


a 5 


h 




x 2 


C3 


d 3 


c 4 


d\ 


C5 


d 5 




X 3 


a 6 


b 6 


a 7 


b 7 


a 8 


b S 




x 4 


C6 


d 6 


C7 


d 7 


C8 


d s . 




_x\ 



where 



f u (x, y) = vo + v\x + v 2 x 2 + v 3 y + v^xy + v^x 2 y + i/ 6 y 2 + v 7 xy 2 + v%x 2 y 2 , 
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for v = a,b, c, d. 

The symmetry condition (iii) reduces the number of unknowns by half since m(x, y) 
becomes 



m(x,y) 



' 1 " 


T 


y 




y 2 




y 3 




y 4 




y 5 





a b a\ bx a 2 b 2 

b s a 8 b 7 a 7 b 6 a 6 

a 3 63 04 64 a 5 65 

65 05 64 a 4 63 a 3 

fl6 bg a 7 b 7 a% bt 

6 2 a 2 61 01 6 «o 





1 " 




a; 




x 2 




X 3 




x 4 




X 5 



For convenience, we denote ^ i/ = a + b. Thus, by (i), we have 



v—a,b 
8 



m 



1(1,1) = 2^] J] 1/4 = 2^(04 + 60 = 1. 
i=0 iy—a,b i—0 

By (ii), we have the following 13 nonlinear equations 



E 

i^— a, 6 

E 

u—a,b 

E 

i^— a, 6 

E 

u—a,b 

E 
E 

i/=a,b 

E 

iv— a, 6 

E 

iv— a. 6 

E 

i/— a, 6 

E 

i/— a, b 

E 

v=a,b 

E 

V — CL. b 

E E**4 

«— u—a,b 



VqVs = 
^2^6 = 

(z/lf6 + ^2^7) = 
{v a V 7 + ViVg) = 
(t/ 2 f 3 + ^5^6) = 
(2^5 + ^3^8) = 
(^0^6 + ^1^7 + ^2^8) = 
{v V 2 + + v 6 v&) = 

(uius + V 2 V/x + V4lS e + U 5 U 7 ) = 
(v l>4 + ^1^5 + ^3^7 + V±V%) = 
(^0^3 + VyV/x + ^2^5 + ^3^6 + ^4^7 + V 5 V 8 ) = 
(U fi +VilS 2 + + + V§V-! + V 7 V$) = 



(1) 



(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
(10) 

(11) 
(12) 

(13) 
(14) 
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The first moment condition (iv) for M — 1 yields the following six linear equations: 

ao + ai + a-2 = bo + h + b 2 (15) 

a.3 + a 4 + a 5 = b 3 + b 4 + b 5 (16) 

a e + a 7 + a$ = b 6 + b 7 + b s (17) 

a Q + a 3 + a 6 = b 2 + b 5 + b 8 (18) 

o,i + a 4 + a-; = bi + b 4 + b 7 (19) 

ai + a 5 + a 8 = b Q + b 3 + b 6 . (20) 



We need to find the otj's and &j's which satisfy the equations (l)-(20) simultaneously. 
We proceed by specifying necessary conditions derived from these equations. 



8 8 

Lemma 2.1. ^]ai = ^ o 4 = -. 



i=0 i=0 

8 8 

Proof. Adding (15)- (17) together, we have ^ = ^ bi. By equation (1), the result 

»=o »=o 

follows. ■ 

Next, equations (2), (3), (8), (9), and (14) imply 

(v\ + (z/i + v 7 f + (V3 + v h f + (i/ + ^2 + ^6 + ^s) 2 ) = (21) 

i/= a, 6 

We use equations (10) and (1 1) to get 

Y (M v o + V2 + V& + v&) + {yi + vi){v3 + vs)) = 0. (22) 

u— a,b 

From equations (4), (5), and (13), we have 

((H + V7)(vo + + + v s ) + 1/4(1* + z/ 5 )) = 0. (23) 

u—a.b 

Also, we have, using equations (6), (7), and (12), 

M^i + ^r) + (vo + V2 + V6 + + = 0. (24) 

u—a.b 

Lemma 2.2. 

1 1 

«0 + a 2 + «4 + «6 + a 8 = 7T H ?= COS a, 

8 4V2 

b + b 2 + b 4 + b 6 + b$ = i H %= sin a. 

8 4^2 



6 



LAI AND ROACH 



Proof. By equation (21) and (22), we have 

(("4 + to + ^2 + ^6 + ^ 8 )) 2 + ((^1 + ^) + to + ^ 5 )) 2 ) = |. 

a, 6 

By Lemma 2. 1 and letting a* = a + a 2 + 04 + o 6 + a 8 and 6* = 60 + &2 + &4 + &6 + &8> 
we have 

or 

(a* - i) 2 + (6* - i) 2 = — . 
v 8 7 v 8 7 32 

Thus, we are able to conclude the proof. ■ 

Let a = a\ + a 4 + 07. Observe that equation (19) implies a = bi + 64 + 07. It follows 
from (21) and (24) that 

V ((1/4 + (H + ^7)) 2 + ((^0 + f 2 + f6 + «*) + to + ^s)) 2 ) = 7J- 
i/= a, 6 

Thus, 2(a) 2 + 2(j — a) 2 = |. That is, we have a = or ^. Similarly, it follows from 
(21) and (23) that 

XI (to + ^4 + ^ 5 ) 2 + (to + f 2 + v% + m) + to + ^)) 2 ) - ^. 

1/— a, 6 

Let a = a 3 + 0S4 + 05 = o 3 + o 4 + 6 5 where the second equality comes from equation (16). 
Thus, 2(a) 2 + 2(± - a) 2 = |. That is, a = or a=\. Therefore, we have the following 
four cases to consider: 

• Case 1: a\ + 04 + a 7 = 1/4 and a 3 + 04 + a 5 = 1/4 

• Case 2: ai + a 4 + a 7 = 1/4 and a 3 + a 4 + a 5 = 

• Case 3: a\ + a 4 + a 7 = and a 3 + 04 + a 5 = 1/4 

• Case 4: ai + a 4 + a 7 = and a 3 + a 4 + a 5 = 

Since Case 3 is a rotation of Case 2, we will only study Cases 1,2, and 4. As will be 
discussed, the different cases are actually conditions associated with the zeros of m(x, y) 
along the axes. 

2.1. Complete Solution of Case 1 

Let us first consider Case 1 where both 01+04 + 07 = i and a 3 + 04 + 05 = i. By 
Lemmas 2.1 and 2.2, we have 

01+03 + 05 + 07 = \ \= cos a, &i + b 3 + b 5 + 67 = ^ \= sin a. 

o 4y2 8 4y2 

However, (ai + o 4 + 07) + (a 3 +04 + 05) = |. It then follows that 

111 31 

204 = 1 = cosa or 04 = 1 = cosa. 

2 8 4^ 16 8^2 
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Similarly, i> 4 = ^ + 7^75 sin a. Let us now summarize the above discussion combined 
with Lemma 2.2 in the following lemma. 

Lemma 2.3. 

3 1 ,31. 

fl4 = l6 + ^7! cosa ' ^le + ^Tl 8111 "' 

a 1+ a r = l--L CO sa, &1 + fe 7 = 1 - ^= sina, 

as + a 5 = ^ - ^ cos a, 6 3 + & B = ^ - ^ sin a, 

ao + a 2 + a 6 + a 8 = --^ H ^cosa, &o + &2 + &6 + &8 = --^ H ^7= sin a. 

16 8V2 16 8V2 

We are now ready to solve for the 18 unknowns. The complete parameterization of Case 
1 is given in Theorem 2.1. 

Theorem 2.1. For any /3,j e [0,2n], let 



a 


= 2(13 


-7)4 


7T 

"4' 


P = 


1 1 

16 8V2 


cos a, 


and q = — \= sin a. 

16 8>/2 


a 


= (- 


-P(l + 


cos(/3 


-7)) 


— q sin(/3 - 


-7)~ 


\/p 2 + q 2 (cos(3 + cos7))/4 




= (" 


-p(l- 


cos(/3 


-7)) 


+ q sin(/3 - 


-7)" 


\Jp 2 + q 2 (cos f3 — cos7))/4 


«6 


= (" 


-p(l- 


cos(/3 


-7)) 


+ q sin(/3 - 


-7) + 


\/p 2 + g 2 (cos (3 — cos7))/4 



a$ = (-p(l + cos(/3 - 7)) - g(sin(/3 - 7) + \/ p 2 + <? 2 (cos/3 + cos7))/4 



60 


= (" 


-q(l + cos(/3- 


-7)) + 


psin(/3 - 


- 7) - vV 


+ g 2 (sin /3 


+ sni7))/4 




= (" 


-g(l — cos(/3 - 


-7))- 


psin(/3 - 


- 7) - Vp 2 


+ g 2 (sin/3 


— sin7))/4 


6 6 


= (" 


-^(l — cos(/3 - 


-7))- 


psin(/3 - 


- 7) + Vp 2 


+ g 2 (sin/3 


— sin7))/4 


6 8 


= (- 


-g(l + cos(/3 - 


- 7)) + psin(/3 - 


-7) + vV 


+ q 2 (sin /3 


+ sin7))/4 


ai 


P 
2 


+ ^"vV + <7 2 


cos/3, 


61 = 




+ q 2 sin /3 




03 


P 
2 


+ ^vV + <7 2 


cos 7, 


h = 




+ g 2 sin 7 




a 5 


P 
2 


- ^"vV + <7 2 


cos 7, 


b 3 = 


9 1 AT 
2 "2^ 


+ q 2 sin 7 




a 7 


P 
2 


- ^Vp 2 + 1 2 


cos/3, 


b 7 = 


2 "5^ 


+ (j 2 sin /3 






1 


-P, h = 


1 










04 


4 


A'* 











then m(x, y) with these coefficients ai and bi satisfies the properties (i), (ii), (Hi), and(iv). 
On the other hand, ifm(x,y) satisfies the properties (i), (ii), (Hi), (iv), and the conditions 
of Case 1, then the coefficients ai and bi can be expressed in the above format. 
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Proof. The proof consists of deriving necessary conditions from various combinations 
of the nonlinear equations which eventually leads to a complete parameterization that 
satisfies each nonlinear equation individually. 

Assume the conditions of Case 1, i.e. a\ + a 4 + a 7 = 1/4 and a 3 + a 4 + a 5 = 1/4. 
Using equations (9) and (14), we have 

J2 ((yo + v if + ("3 + + K + ^s) 2 + v\ + + vj) = 

u— a,b 

Using (13), the above equation becomes 

^2 {(vo + v\ + v 2 f + (v 3 + + v 5 ) 2 + (i/ 6 + v 7 + z/ 8 ) 2 ) = i. 

a, 6 

Since 03 + a 4 + a$ = 63 + 64 + 65 = ^, we have 

aa + ffli + a 2 = 0, a 6 + a 7 + a 8 = 0, 60 + ^1+^2 = 0, 6 6 + b 7 + b s = 0. (25) 
Similarly, by equations (8) and (14), we have 

V (fa + ^ 6 ) 2 + + v 7 f + (z/ 2 + v 8 f + 4 + + vl) = \. 
— \ 8 

v— a,b 

Using (12), the above equation becomes 

^ ((^0 + ^3 + ^e) 2 + {vi + v 7 f + {yi + v b + is 8 f) = i. 

u—a : b 

Since a\ + a 4 + a 7 = &i + 64 + fe 7 = ^, we have 

ao + a 3 + a 6 = 0, a 2 + a 5 + a 8 = 0, 6 + h + h = 0, 6 2 + &5 + &s = 0. (26) 

We now use (25), (26), and the equations in Lemma 2.3 to simplify the equations (2)-(14). 
By (25), equation (13) becomes 

= + "ate + - ^) 

= + v ^ 2 + ~ ^ + 1/11/7 

b>—a,b 

u—a,b 
u—a,b 

since J2v= a b(~ (i — ^ 4 ) 2 + ^(i — ^ 4 )) = by Lemma 2.3. A similar situation exists 
for equation (12), thus equations (12) and (13) simplify to 



0305 + hh = 0, aia 7 + bib 7 = 0. 



(27) 
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We are now able to solve for a 1 ,a 3 , a 5 , a-?, b\, 63, 65, and 67. For simplicity, we denote 
a\ + 0,7 — p and b\ + bi = q. By (27) we have 

a 2 -pa 1 + b 2 -b 1 q = or (a, - |) 2 + (h - |) 2 = \{p 2 + q 2 ). 
Thus, we have 

V , 1 



ai = 2 + 2^ +?2 C0S/3 ' &1 = 2 + 2^ p2 + q2 Sin ^' 
a 7 = I - TjVp 2 +Q 2 cos/3, fe 7 = I - t^vV + <7 2 sin^. 

Similarly, we have 

a 3 = | + ^vV + <? 2 cos 7 , 63 = I + TjVp 2 + q 2 sin7, 
a 5 = | - t^vV + <? 2 cos 7 , 6 5 = I - ^vV + q 2 sin7- 

Next we simplify the equations (2)-(l 1). Note that the addition of equation (10) and (1 1) 
is 

<>4(>o + V2 + v§ + vs) + (vi + vr){v3 + ^5)) = 0. (28) 

v— a,b 

After substitution of the equations in Lemma 2.3, the left hand side of equation (28) 
becomes 

2(J -P)(-P) + 2(-p)(-p) + 2(1 - </)(-«) + 2(-g)(-«) 
= 4[p 2 - -p + (— ) 2 + g 2 - -q + (— ) 2 ] - 8(— ) 2 

= 4[(p-— ) 2 + (<?-— ) 2 ]- — 
16 ; v 16 ; J 32 

111 

= 4[— — - cos 2 a H — — sin 2 a] — — = 0. 
L 128 128 J 32 

That is, the equations (10) and (11) are linearly dependent and consequently only one 
needs to be considered. We will deal with equation (10) later. Furthermore, for equations 
(4)-(9), we have the following 

{viVg + v 2 v 7 ) = ^2 vivs + v\v-! + v 2 v§ + v 2 v 7 

v=a,b u—a.b 

= ^2( l/ i + v 2)(y% + vi) = ( _t, o)(-^s) = v ° v& > 

u—a.b v=a,b u—a,b 



^ [v^vi + v\v%) = Y, 1/01/7 + 1/01/8 + 1/11/7 + 1/11/8 



v—a.b u—a.b 



v—a,b v—a.b v—a.b 
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v—a.b v—a.b 

= ^2 {"2 + V5){V3 + V 6 ) = ^2 {-V8){-Vo) = ^2 VqV 6 , 
v—a.b i/=a,b v=a,b 

^ {VQVs + V3I/&) = ^0^5 + + ^3^5 + ^3^8 

v—a,b v—a.b 

= J! ("o + ^3)(^5 + ^s) = X! (-^eX-i^) = X! I/2l/6 ' 

v—a.b v—a,b v—a,b 

^ (^0^6 + ^1^7 + ^2^8) = ^0^6 + ^2^8 = ^0^6 + ^0^8 + ^2^6 + ^2^8 

v—a.b v—a.b 

^2 ( v o + V2)(y % + v s )= ^2 (-^iX -1 ^) = ^2 l,lL " 7 ' 



v—a.b v—a.b v—a.b 



v—a,b v—a.b v—a.b 



^ {y^2 + ^3^5 + ^6^8) = X] 7/07/2 + 7/07/8 + 7/27/6 + 7/67/8 
v—a,b v—a.b 

= ^2 {"0 + v%)(y 2 +v 6 )= ^2 (- jy 3)(-^5) = 

u—a.b v=a,b v—a,b 

Therefore, only $^„_ a 6 ^o^s = and X^= a 6 ^2^6 = remain to be solved in addition to 
(10) since equation (14) follows from (1)-(13). 

In order to solve equations (2) and (3), note that (25) and (26) imply a = a 5 + a 8 — 01 
and b = b 5 + b$ — b\ . Thus, (2) becomes 

a\ + a 8 (a 5 - ai) + b 2 s + b s (b 5 - h) = (29) 

while (3) becomes 

(a 5 + a 8 )(a 7 + o 8 ) + (65 + &s)(&7 + b 8 ) = 

or 

ag + a 8 (a 5 + a 7 ) + bj + b$(b 5 + b 7 ) = -a 5 a 7 - 6567. (30) 
The subtraction of (29) from (30) yields 

a 8 (ai + a 7 ) + b s {h + b 7 ) = -a 5 a 7 - b 5 b 7 . 

So, we have 

pa 8 + qb s = -a 5 a 7 - b b b 7 . (31) 

The addition of (29) and (30) yields 

2 / a 7 ~ a i \ 7 2 /j b 7 — bi 1 iin 

a 8 + («5 H 2 ) as + "s + (°5 H 2 — > 8 = _ 2^ a5<X7 + ^ 
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which is 

(as + \{a 5 + ^p)) 2 + (6» + 1(6 B + ^i)) 2 = R, (32) 
for some known value R\. In addition, equation (31) can be rewritten as 

P(a 8 + \{a 5 + ^=^)) + q(bs + \{h + ^y^)) = R* (33) 

for some known value R 2 . Equations (32) and (33) can be solved simultaneously. Thus, 
we obtain the expression for the ctj's and bi's given in Theorem 2.1. 

Finally, to satisfy (10), we put these a, and &j into (10) and simplify the equation yielding 

cos(/3 - 7) = -j= cos(a - j3 + 7) + sin(a - /? + 7). 

Solving, we find that a = 2(/3 — 7) + 7r/4. 

The above discussion shows that if m(a;, y) satisfies the properties (i), (ii), (iii), (iv), and 
the conditions of Case 1 , then its coefficients a, and &, can be expressed as the two-parameter 
family given in the statement of Theorem 2.1. 

On the other hand, to verify that m (x, y) with the coefficients a^and bi satisfies the proper- 
ties (i), (ii), (iii), (iv), and the conditions of Case 1, we just substitute the solutions back into 

equations (l)-(20). This completes the proof. ■ 



2.2. Solution of Case 2 

In this section, we consider Case 2 where a\ + + a 7 = 1/4 and a 3 + 04 + a 5 = 0. 
By Lemmas 2.1 and 2.2, we have ai + a 3 + a 5 + a 7 = | — cos a. So, we have 

11 1 1 

a 4 = 1 = cos a and 64 = 1 = sin a. 

16 8^ 16 8V2 

By (9),(13), and (14), we have 

V] ((i>o + v 2 ) 2 + (^3 + ^4 + ^) 2 + (^6 + ^7 + ^s) 2 ) = 77- 

v— a,b 

The assumptions of Case 2 imply 

(Of) + V\ + V 2 f + K + ^7 + ^s) 2 ) = 

i/— a, 6 

i.e. 

(a + 01 + a 2 ) 2 + (ae + «7 + a 8 ) 2 + (60 + &i + ^) 2 + (&b + &7 + &s) 2 = 77- 

o 

By (15) and (17), we have 

(o + ai + a 2 ) 2 + (a 6 + a 7 + a 8 ) 2 = 
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By Lemma 2.1, the above equation becomes 

I i 

(4 - K + a 7 + a 8 )) 2 + (a 6 + a 7 + a 8 ) 2 = — . 

It follows that a 6 + 07 + o 8 = or a 6 + 07 + a 8 = j. Thus, Case 2 branches out into two 
subcases. 

1 1 

• Subcase 2a: 01+04+07 = -, 03+04 + 05 = 0, 06+a 7 + a 8 = 0, a + ai+a 2 = -. 

• Subcase 2b: 01+04+07 = i, 03+04+05 = 0, a 6 +a7+a 8 = ^, ao+ai+a 2 = 0. 

We only consider subcase 2a here. The subcase 2b can be treated similarly and is left to 
the interested reader. Theorem 2.2 gives the complete solution of Subcase 2a. 

Theorem 2.2. 

1 ) For any 7 <G [0, 27r], 03 = 63 = a 4 = b 4 = a 5 = 65 = 0, 

3 1 11 

ai = — -= cos 7, 07 = — H -= cos 7, 

16 8V2 ' 16 8V2 

o = ^1 + V2cos7 ± + V2(cos7 + sin 7)^ , 
77^ ^1 + \/2cos7 + ^2 + V2(cos 7 + sin 7)^ , 



a 2 = 



a 6 



= -7- ( -1 - V2COS7 + \Z 2 + V2(cos7 + sin 7) ) , 



32 



a 8 = ^-1 - \/2cos7 ± y^2 + V2(cos7 + sin 7)^ , 

31 1 1 
61 = — =sin7, 07 = — H =sm7, 

16 8\/2 ' 16 8V2 



77^ ^1 + \/2sin7 + ^2 + \/2~(cos 7 + sin 7)^ , 
6 2 = 7^- ( 1 + \/2shi7± v/2 + \/2(cos7 + sin 7) ) , 



^6 = ^-1 - V2sin7± ^2 + V2(cos7 + sin 7)^ , 
6 8 = 77^ ^-1 - V2smj + ^2 + V2(cos7 + sin 7)^ . 

2) For any a G [0, f ) U (f,27r], Zef 



7 = 



-ia + ^ or ±a if 0<a<f 

|a+| or -§a-f if f < a < 2tt 

11 11 

p = -— — cos a, q = — — -= sin a, 

16 8\/2 16 8\/2 
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a = --7- (-1 + 8p + s H — [(-s - 8p + 8q) cos-f + ss'm-f] 

62 \ p — q 

02 = 777: ( 1 — 8» + sH [(— s + 8p — 8q) cos 7 + ssin7l 

32 \ p — q ' 

a 6 = --7- (l + 8p + s-\ — Us + 8p- 8q) COS7 - ssin7] ] 

32 V p-q J 



<k = -7; ( — 1 — 8p + s H — [(s — 8p + 8g) COS7 — sshi7] 



32 



b = -7- ( 1 - 8q + s H — [(s + 8p - 8g) sin7 - s cos 7] 

62 \ p — q 

b 2 = -7-- ( -1 + 8q + s H — [(s - 8p + 8q) sin7 - s cos 7] 

32 v p-q 

K = f -1 - 8g + s H — [(-s - 8p + 8q) sin 7 + s cos 7] 

32 V. p-q 

b s = --7- M + 8<? + s H - — [(-s + 8p - 8q) sin7 + SCOS7] j 



fll = 


— (3 + 8p-8icos 7 ), 
Id 


h = 


4(3 + 8<z-8tsin7) 
Id 


a 7 = 


^(l + 8p + 8icos7), 


b 7 = 


— (1 + 89 + 8*8^7) 


«3 = 




b 3 = 


^q-s), 


a 5 = 


^(8P-), 


b 5 = 


±(Sq + s) 


C14 = 


~P, 


b 4 = 


-q- 



3) For a 








1 


1 


1 


1 







1 


1 


2 


2 


1 


1 
























-1 


1 








1 


-1 


1 





-1 


1 


1 


-1 





1 





-2 


2 


2 


-2 



















or — 














8 





-1 


1 


1 


-1 





16 





-2 


2 


2 


-2 



























-1 


1 








1 


-1 







1 


1 


1 


1 







1 


1 


2 


2 


1 


1 



If m(x, y) has the coefficients given by 1), 2), or 3), then m(x, y) satisfies the properties 
of (i), (ii), (Hi), and (iv). Conversely, if m(x, y) satisfies the properties of (i), (ii), (Hi), (iv), 
and the conditions of Subcase 2a, then the coefficients of m(x, y) can be expressed as 1), 
2), or 3). 

Proof. Assume the conditions of Subcase 2a, i.e. 

ai + a 4 + a 7 = a 3 + a 4 + a 5 = 0, a 6 + a 7 + a 8 = 0, a + ai + a 2 = -V 
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We immediately have from Lemma 2.2 



11 ,11 

a4 = — H = cos a, 64 = — H = sin a. 

16 8V2 16 8^2 



By equations (8), (12) and (14), we have 

^ (>o + ^3 + ^e) 2 + {vi + v-if + (v> 2 + + vg) 2 = \. 

v—a,b 

It follows that 

ao + «3 + «6 = 0, 6 + 63 + be = 0, 

a 2 + a 5 + as = 0, 6 2 + b 5 + b 8 = 0. (34) 

We now simplify the 13 nonlinear equations (2) -(14). With (2), (3), and (34), we may 
simplify (9) as follows: 

= {vqV2 + ^3^5 + V&V%) = ^ (^0 + ^6)(^2 + ^s) + ^3^5 

z/— a, 6 v—a.b 

= ^2 ( _ ^)(-^5) + ^3^5 = 2 J! ^ 35 ^ 
v—a.b v—a.b 

Recall a 3 + a 4 + a 5 =0 and 6 3 + o 4 + o 5 = 0. We can now solve for 03, as, 63, 65. 
Indeed, (35) can be rewritten as (04 + a 5 )a 5 + (64 + 65)65 = 0. After simplifying, we have 

"4.2 °4\2 _ 1/ 2 j2\ 



Thus, we have 



03 = 2 + 2 vV + ^cos/?, 63 = - + -y/p 2 + q 2 sin/3 
a 5 = I - \ Vp 2 + <7 2 cos /3, & 5 = I - ^vV + <? 2 sin/3, 
where p = -a 4 and q = —64 (different p and q from section 2.1). We note that 

p 2 +q 2 + \(p + q) = 0. (36) 
Similarly, we may simplify (8) to be 

= ^ {v\v 7 + v q vq + v 2 v%) = v ii> 7 + (v a + v 2 ){vti + v?,) 

v—a,b v—a.b 

= 1/11/7 + (4 ~ ^iX-^) = 2 H (^^7 - g^r)- 

v—a.b u—a,b 



(37) 



Recall ai + a7 = i — a 4 and 61 + 67 = ^ — 64. We may solve for a\ and b\ to get 
ai = I — 04 — 07 and 61 = | — 64 — 67. Putting them into (37), we have 

1 \ a 7 ( 1 , , \ , 67 „ 

— — 04 — a-j \ a-j — — + — — 04 — 67 67 — — = 



4 7 8 V 4 
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or 



a 7 - 



1 - 2 

a 4 



+ [b 7 - 



i-h\ 2 _(a 4 -i) 2 + (& 4 -iy 



It follows that 



a 7 = 


? + 


l 

16 + 


b 7 = 


h 


1 

16 + 


ai = 




3 

16 ~ 


bi = 


h 


3 

16 ~ 



8 7 \ q 8 ' ( '° S " ' 



sin 7, 



We add the left-hand side of (10) and (11) together to get 



= 2 ^ 4 

y= a, 6 

= 2 1 2a? - t< 



2[v 1 + v 7 



z/— a, 6 



|- -a 4 + 2 26; 



0. 



That is, only one of (10) and (11) needs to be solved. Thus, we deal with (10) later. 
Turning our attention to equation (13). We have by (35) and (36) 

u—a.b 

= ^2 ("o + v-^j V\ + ^1/3 + IS 5 ^j V A + V 7 ^Vq + IS 8 ^J 
= ^2 -2v x v 7 - v\ - v 7 - v\ + j {Vx + V 7 ) 



v—a.b 



E (-(^1 + ^7) -^f + ^1+^7)) 



y—a,b 
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- E (i-^) 



2 1/1 

2 1/1 



vA = 0. 



x 4 "V ' " 4 4V4 

i/—a,b 

That is, (13) holds for these 04, 0,1,0,7, b\, 67. Similarly, (12) is satisfied in that 

^Z ( ^0^3 + ^1^4 + ^2^5 + ^3^6 + ^4^7 + ^5^8 ) 

v=a,b ^ ' 



E ( -^3+^(1-^4) -^ 2 ) 

'=0,6 ^ \ / / 



u—a,b 



= -E(^-t)=°- 

z^— a, 6 



We now show that equations (3), (9), (34) and (35) imply (6). The addition of (6) and (9) 
yields 

^ ^6 ( ^5 + ^8 ) + »2 ( ^0 + ^3 ) + ^3^5 = ~ 2 " 2 ^ + ^5 = 0. 

v=a,b V / V / j/=a,b 

Similarly, the equations (2), (9) and (35) imply (7). Since (8) holds for these ai,br,ar, by 
under the assumptions of (2) and (3), we further simplify (4) by adding (4) and (8) together. 
That is, 

= ^2 V 2 V-! + V X V(, + VyV-j + V Vq + ^2^8 
v— a,b 

^2 ^ 6 ( U ° + Vl ) + V2 ( Vl + 1/8 ) + Vl V7 

v=a,b ^ ' ' 



v—a,b 

= E (7^6+^7)= E ( i^ + ^7 ). 



h x 4 V ^A4" 8 

v=a,b x ' v—a,b 



Similarly, the sum of equations (5) and (8) is equivalent to 



E ("i^ + \ v *) = E + \ v *) = °- 

v=a,b ^ ' v—a,b 



(38) 



However, equations (37) and (38) are equivalent by using a,Q+aj+a% = and b^+bj+bg, 
0. Thus, we only need to consider one of them. 

In summary, we only need to solve the following equations 



aoa 8 + b Q b 8 = 0, a 2 a 6 + b 2 b$ = 0, 
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^ {V7 + 2^8) =0, ^ (^0^4 + ^4^8 + ^1^5 + V$Vt) = 0. 
v—a.b u—a,b 

Using the linear relationships, we have 

a = i-ai+a 5 + a 8 , 6 = ^ - 6 X + b 5 + b 8 

a-2 = -05 -ag, 6 2 = -65 - &8, «6 = -«7 - as, h = -b 7 -b s - 

Putting these linear equations in (2) and (3), we get 

a\ + (a 5 - ai + ^)as + &s + ( b 5 ~h + = 

a\ + (a 5 + a 7 )a 8 + a 5 a 7 + &§ + (h + M & 8 + &5&7 = 0. 

Subtracting the first one of the above two equations from the second one, and using 

ai + a 7 = p + |, &i + 6 7 = q + \, we get 

pa 8 + g& 8 = -a 5 a 7 - b 5 b 7 . 
Using these linear relationships and (1 1), we get 

pa$ + qb 8 = i ^ (i/!i/ 5 + 1/31/7(1/5 - ^1 + | W)- 

a,b 

It follows that 

(21/51/7 + ^1^5 + ^3^7 + (^5 + = 

1/— a,b 

which can be simplified to be 

E (7^ + ^) = 



'4 

z^— a, 6 



That is, we have 

= 04 + 64 + ^(05+65) 

= -(a 4 + 64) + -(p + o + -vV + q 2 (cos/3 + sin f3) 
o 4 z 

= -Lv?+^sin(J+/?). 

It follows that either [3 = -\ , ox p 2 + q 2 = which the latter occurs when a = -^f. 

We first consider a = — ^f. In this situation, we have p = and q = 0. It follows that 
03 = 63 = a 4 = 64 = 05 = 65 = 0. It is clear that pa 8 + qb$ = — J2u=a b v ^ holds. 
Thus, we only have two equations to solve 

a a$ + b b 8 = 0, a 8 + b$ = ~(d7 + h)- 
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Using a — j — ai + a 8 and b a = \ — b\ + 6 8 , we solve the above equations and get 

a 8 = ^ — 1 — V2 cos 7 ± y^2 + \/2(cos 7 + sin 7)^ 



h 32 



-^^-l-\/2sin7=F y 2 + \/2(cos 7 + sin 7) 
a ° = ^2 ( 1 + ^2cos7 ± \jl + \/2(cos7 + sin7)^ 
b = ^1 + \/2sin7 =F \/ 2 + \/2(cos7 + sin 7)^ . 
Using the linear relationships, we have 

a 2 = ^1 + \/2cos7 =F y^2 + \/2(cos7 + sin 7)^ 
^^1 + ^2 sin 7±y / 2 + \/2(cos 7 + sin 7)^ 



62 = 



a 6 



— [ - 1 - \/2 cos 7 =F v2 + V2(cos 7 + sin 7 



32 



6 6 = ^ — 1 — \/2 sin 7 ± \J 2 + \/2(cos 7 + sin 7)^ 

3 cos 7 3 sin 7 1 cos 7 1 sin 7 

fll = l6"^7f' 6l = l6"^7f' a7 = l6 + ^7f' 67 = l6 + ^7f' 

Next we consider (3 = — f or /3 = We only need to consider /3 = — j while 
p 2 + g 2 7^ because the other case is a rotation of this one. For /3 = — j, we have three 
equations to solve: 

0008 + 60^8 = 0, pa$ + qb$ = - ^2 a 5 a 7 , a 8 + b 8 = -^(a 7 + b 7 ). (39) 

Assuming p 7^ q, i.e. a ^ |, we can solve the second and third ones in the above three 
equations for a 8 and b$: The solutions for a 8 and 6 8 are 

6 8 = — !— ( - ^(a 7 + b 7 ) + V a 5 a 7 

a 8 = -\{a 7 + 67) + — !— ( - ^(a 7 + 6 7 ) + a 5 a 7 ) 
^ p g\ ^ u=ab J 

leaving the relationship between a and 7 as a a 8 + 6 °8 = 0. Upon substitution and 
simplification, we have 



sin(o; + ^) - 2sin(7+ ^)y2(l - sin(a + ^) + sin(2 7 ) - 2 = 0. (40) 

This equation (40) has a one parameter family of solutions given by 

_{-\a+ 7 -f or \a-^- if 0<a<f 



CONSTRUCTION OF BIVARIATE SYMMETRIC WAVELETS 



19 



The final case is when j3 = — \ and p = q (i.e. a = f) which implies that 

on = 61 = a 4 = hi = \, a 3 = a 7 = b 5 = b 7 = 0, a 5 = 63 = ~- 

5 o 

This reduces (39) to 

a a 8 + b b s = 0, a 8 + b s = 0, 

which has two solutions: a 8 = b 8 = or a 8 = and fe 8 = — jg yielding two rational 
solutions: 








1 — 1 


1 


1 


I— 1 







1 


1 


2 


2 


1 


1 
























-1 


1 








1 


-1 


1 





-1 


1 


1 


-1 





1 





-2 


2 


2 


-2 





8 





-1 


1 


1 


-1 





' 16 





-2 


2 


2 


-2 



























-1 


1 








1 


-1 







1 


1 


1 


1 







1 


1 


2 


2 


1 


1 



2.3. Solution of Case 4 

Finally, we consider Case 4 where both ai + a 4 + a 7 = and a 3 + a 4 + a 5 = 0. From 
equations (16) and (19) we have 6 3 + & 4 + &5 = and fei + & 4 + 67 = 0. As before, using 
equations (9), (13) and (14), we have 

V {{ya + v\ + V2? + (^3 + "4 + Vbf + Oe + V7 + y&) 2 ) = \- 

u ° 
v=a,o 

It follows that Ylv=a,b(( l 'o + + v 2 ) 2 + {vq +V7 + v&) 2 ) = |. Using equations (15) and 
(17), we have 

(o + 01 + a 2 ) 2 + (a 6 + a 7 + a 8 f = 
By Lemma 2.1, the above equation is 

(i - (o 6 + a 7 + a 8 )) 2 + (a 6 + a 7 + a 8 ) 2 = 

Thus, + S7 + a 8 = i or a 6 + a 7 + a 8 = 0. 
Also, using equations (8), (12) and (14), we have 

((^0 + v z + i^ 6 ) 2 + (i/i + 1/4 + z/ 7 ) 2 + (i/ 2 + ^5 + ^s) 2 ) = ^. 

a,o 

It follows that X)i/= ,6 (("0 + v 3 + v§) 2 + (i>2 + V5 + ^s) 2 ) = |. So, we have 

(a + a 3 + a 6 ) 2 + (a 2 + a 5 + a 8 ) 2 = 

lb 

Again by Lemma 2.1, (| - (a + a 3 + a 6 )) 2 + (a + a 3 + a 6 ) 2 = ^. Thus, a + a 3 + a 6 = 



or a + a 3 + a 6 = \. 
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Therefore, we have four subcases to consider. In addition to 

01 + a 4 + a 7 = 0, a 3 + a 4 + a 5 = 0, 6i + 6 4 + 67 = 0, 63 + 64 + 65 = 0, 

we have 

• Subcase 4a: 

Oo+ai+a 2 = i, ao + o 3 + a 6 = i, a 6 + a 7 + a 8 = 0, 02 + 05+03 = 
6 + 6i + 6 2 = i, 6 + &3 + &6=0, 6 6 + 6 7 + 6 8 = 0, & 2 + 6 5 + 6 8 = l 

• Subcase 4b: 

1 1 

a + ai + a 2 = 0, a + a 3 + a 6 = -, a 6 + a 7 + a 8 = -, a 2 + a 5 + a 8 = 0, 

6 + 6i + 6 2 = 0, 6 + 6 3 + fo 6 =0, 6 6 + 6 7 + 6 8 = ^, b 2 + h + b 8 = ^ 

• Subcase 4c: 

00+01+02 = ^, ao + 03 + 06 = 0, a 6 + a 7 + a 8 =0, a 2 + a 5 + a 8 = ^ 

&o + &i + &2 = 7, 60 + 63 + ^6 = 7, 65 + 67 + 68 = 0, 62 + 65 + 68 = 0, 



4' 



• Subcase 4d: 



ao + 01 + a 2 = 0, a + o 3 + a 6 = 0, a 6 + a 7 + a 8 = i, a 2 + a 5 + a 8 = ^ 
6o + 6i + 6 2 = 0, 6 + 6 3 + & 6 = i, 6 6 + 6 7 + 6 8 = i, & 2 + 65 + b 8 = 0. 

We only study the Subcase 4a and leave the other three subcases to the interested reader. 
With the linear constraints, we tackle the nonlinear conditions (2)-(14). We use (2) and (3) 
to simplify (8) and (9) as follows: 

= ^ + v\v 7 + v 2 f& = (fQ + ^>)(+6 + ^s) + ^1^7 

u—a.b u—a.b 



(41) 



v—a.b x 7 u—a,b 

and similarly, 

= ^2 (^2 + ^5 + v d v &) = ("° + + ^s) + ^3^5 

u—a.b u—a.b 

= o 3 )(-o 5 ) + o 3 o 5 + (-6 3 )(i - 65) + 6365^ (42) 

— -^+0305 + 6365). (43) 
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As we did before, we have a 4 = — yg + cos a and 6 4 = — yg + ^= sin a. It follows 
that 

<*4 + 64 + + = 0. (44) 
8 

This fact will be used later. 

With (8), i.e., (41), we can see that (12) holds. Indeed, the left-hand side of (12) is, by 
using (44), 



^2 (^("o + v i) + ^4(^3 + V5) + vi{v§ + v 8 )) = y~] {vijj ~ v i) -aj- a 7) 

'—a,b v—a,b 

= E (\(~ U4 - v t) -al-aj-af) = - ^ (-^4 + {v\ + vif + a\) 

u—a,b i/=a,b 

= - E (j^ + 2v ^ = °- 



'4 

u—a.b 

Similarly, we can show that with (9), i.e. (43), equation (13) holds. 
If we add equations (10) and (11) together, we have by (44) 



^ ^4(^0 + v 2 ) + V x (v 3 + V 5 ) + ^7(^3 + V 5 ) + V±(vq + v 8 ) 
u—a,b 

= ^2 (^(^ - Vl) - V\V± - "7 "4 - V^V-i) 
y—a,b 

^ - 2^4(^1 + 1/7) = ^ jiy 4 + 2v\ = 0. 



4 -v - ^4 

v—a : b u—a,b 

That is, the addition of (10) and (1 1) is always true. We only need to consider one of these 
two equations. Next we simplify equations (4) - (7). 
Adding (4) and (8) with (3), we have 

= (^2(^7 + ^s) + (fp + v\)v<> + V1V7) 

v—a,b 

= E (^+^7) = E (r 6 + i)- 

a,b v—a.b 

Adding (5) and (8) together, we have 

= {V\ + ^ 2 )^8 + V{){V7 + ^6) + V\V7 



(45) 



v—a.b 



^(-1/8+^1/7) = y^ (i^ 8 + iz/ 8 



v—a,b v—a,b 



(46) 



It is easy to see that equations (45) and (46) are equivalent by using a 6 + a-j + a 8 = and 
&o + &i + b 2 = 0. 

Adding (6) and (9) together, we have 



° = E (^3(^5 + ^s) + ^2(^0 + ^3) + 



u—a,b 
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11 11 

= -a 2 + a 3 a 5 + -b 6 + 6365 = -{a 2 + 6 6 ) + -(a 5 + b 3 ). (47) 
4 4 4 o 

Adding (7) and (9) together with (2), we have 

= ^ (^8(^3 + Vq) + ^0(^2 + ^5) + ^3^5) 
v— a,b 

= -a 8 + a 3 <z 5 + 6365 + -6 = 7(08 + 60) + ^(a 5 + 63). (48) 
4 4 4 o 

In fact, equations (47) and (48) are equivalent by using a 2 +a 5 +a 8 = and 60+63+66 = 0. 

Next we solve (41) using ax + a 4 + a 7 = and 61 + 64 + 67 = 0. Letting p = — a 4 and 
q = -64, we have pa 7 - a 2 7 - ^ + qb 7 - tf - \ = 0. That is 

Recalling (44), i.e., p 2 — \p + q 2 — \q = 0, we have 

2+ 2 = l { l- p - q)=1 8{l + ^ {coaa+s[na) 



= ei ( 1 + 7i (cosa + sina) ) = ^( 1 + sin(a + i ) )- 



Thus, we have 



a7= 2"l6 + 16V + Sm(a + 4 } C ° S/3 ' 67 = 2 ~ 16 + 16 V 1 + Sm(a + 4 } Sm/3 ' 
Similarly we can solve (43) using a 3 + a 4 + a 5 = and 6 3 + 6 4 + 65 = giving us 



05 = \ h + h v 1 + sin(a + 1 } C0S7 ' 63 = I - ^ + ^ v 1 + sin(a + i } sin7 - 

It follows that 

01 =f- 1 ^ + nv /1 + sin ( a + f ) cos & 6 i = § ~ 11 + iW 1 + sin ( a + f) sin P 
fl 3 = 2 + Te + TeV 1 + sin(a+ f )cos 7 , &5 = § + Tg + TgV 1 + sin ( a + f) sin ^- 

So, we still need to satisfy (2), (3), (46), (48), and (1 1). Using the linear relationships for 
Subcase 4a, (2) and (3) become 

4 + (J - ai + a 5 )a 8 + 6|(-6i + 6 5 )6 8 = (49) 
ag + (a 5 + a 7 )a 8 + 6g + (65 + 67 - i)6 8 = -a 5 a 7 - 6 5 6 7 + ^. (50) 

After subtracting (49) from (50) and using a x + a 7 = p and 61 + 67 = q, we can replace 
one of these equations with 

(p- ^)a 8 + (q- ^)6 8 = \b 7 - a 5 a 7 - 6567. (51) 
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Moreover, the linear relationships combined with (11) yield 



+ qb 8 = ^a 4 + ^ ^ ^1^5 + ^3^7 + {"5 - vi)v±. 



pas ..... 2 

v— a,b 



Thus, (46), (48), and (11) can be replaced by 



a 8 + 6 8 = -l(a 7 + b 7 ) (52) 
a 8 + 6 8 = -I( a5 -2&i+ 6 3 + 26 5 ) (53) 
a 8 + 6 8 = -67 + ia 4 + 2 ^2 21/51/7 + ^1^5 + + (^5 - ^1)^4- (54) 



Now, we only need to satisfy (49), (51), and (52)-(54). 
Equating the right-hand sides of (52) and (53) gives 



7T 

1 + sin(a + — )(cos7 — sin7 — cos j3 + sin j3) = 0. 

This constraint is satisfied whenever a = — ^p, 7 = /?, or 7 = —j3 — -|. 
Equating the right hand sides of (52) and (54) gives 

~(a 7 + 67) = -67 + ia 4 + 2 ^2 (^Oi + vt) + O3 + ^5)^7 + - ^1^4 

= -67 + 75 a 4 - 2 X! + ^ 



1/— a, 6 



Thus, 



1 9 ^4 v 7 



u—a,b 



1 / 7T 

= 7^Jl + sin(a + -)(cos^-sin^). 

Therefore, we only need to solve (49), (51), and (52) when a = -^,/3=f,or/3 = -^ 
and 7 = f3 or 7 = — f3 — |. 
We begin with a = — then 



a 4 = 64 = -i a 5 = a 7 = 6 3 = 67 = 0, a! = a 3 = 61 = 65 = \ 
o o 

which reduces (49), (51), and (52) to 

(a 8 + i)a 8 + 6 8 = 0, a 8 + 6 8 = 0. 

o 

These equations yield two solutions: a 8 = 6 8 = or a 8 = — and 6 8 = 
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Now, with (3 = 7 = j, we solve for a 8 and b 8 using the linear equations (51) and (52), 
we have 

«8 = ~\ - + sin(a + tt/4)), &s = -| + 

Plugging these solutions into (49) and simplifying we have 



7T / 7T 

cos(a — — ) — 1 = 4W2 + 2cos(a - — ) 



which has no real solution. 

Similarly for = 7 = -?f, (51) and (52) yield 

V2 



«s = -| + — v 7 ! + sin(a + tt/4)), 6 8 = -- + — , 



but now (49) reduces to 



cos(a - ^) - 1 = -4W2 + 2cos(a - ^) 



which has two solutions: a = j ± cos 1 (17 — 8\/5). 

Finally for /3 = | and 7 = — the linear equations produce 



as = — 4 , b 8 = — 4 + --—Jl + S m^-^, 



and similarly for /3 = — ^ and 7 = f 



as = -4, & 8 = -4 + ^ + ^-V 1 + Sin(a "^ 



Both of these choices for a 8 and 6 8 after being substituted into (49) require 



p 2 + q 2 = 



which is satisfied by a = ^. 

Therefore, the complete solution for Subcase 4a has 6 solitary solutions. The four 
rational solutions are given below: 



1 



1 
1 



1 

1 

1 





1 
1 



1 



-1 
-1 



1 
1 

-1 




-1 
1 



1 1 



-1 1 

-1 1 



1 1 

1 1 

-1 1 




-1 1 

1 1 



1 

16 



1 1 

1 -1 

2 
2 
1 -1 
1 1 



2 


-2 
-2 

2 



1 1 



-1 



-1 
2 1 



1 1 

-1 1 



1 


-10 
-10 0-11 



10 11 
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3. ORTHOGONALITY 

In this section, we discuss the orthogonality of the solutions from Cases 1-4. We 
begin with a review of the Lawton condition as well as a well known necessary and 
sufficient condition for orthogonality. We conclude this section with a numerical experiment 
consisting of a one-level decomposition of two gray-scale images using a nonseparable filter 
from Subcase 4a, and compare its performance with Haar and D4. 

Let 



m(e i0Jl ,e i0J2 ) =J2 



'^k,l<^ r 
k,l 

where the hk/s are the a/s and 6j's as discussed previously in Section 2. Define 



(j>(uii,ui2) = ]^[ to (e~, . (55) 



fe=i 

For the coefficients as defined in Section 2, we know that <p is we U defined and is in 
L 2 (R 2 )- Define 

at lt t 2 = / (p{x,y)(j)(x - £i,y - £ 2 )dxdy. (56) 
Jr. 2 

Thus, if ag 1 j 2 = 5g 1 j 2 , then is orthonormal. By the refinement equation 

5 

(j>(x,y)=4 h klM (j>{2x-k l ,2y-k2). (57) 

fci,fc 2 =0 

Using (57) in (56) we have 

O^iA = 4 I hk 1 ,k 2 hk 1 +n 1 -2e 1 ,k 2 +n2-2e 2 I a n u n 2 - (58) 

"i,>i2 \ki,k 2 ] 

Because the supp(cf)) C [0, 5] 2 , the only possible nonzero ae lt t 2 are for — 4 < £i, £ 2 < 4. 
Let a be the vector of length 81 consisting of the a£ 1 .£ 2 , s for some fixed ordering of the 
indices in the range of — 4 < l\, l 2 < 4 and define the matrix 

^(«iA),(ni,n 2 ) — 4 X] hk 1 ,k 2 hk 1 +n 1 -2e 1 ,k 2 +n 2 -2e 2 (59) 
kl,k 2 

for this same ordering. Then equation (58) says that a is an eigenvector of A with 
eigenvalueA = 1, i.e., a = Act. Now, condition (i) of Section 1 implies that 

4 ^ hktMhkt^hM-^hX^y 232 = 1, 

k\,k 2 

i.e. 



4 ^ h ku k 2 h kl -2j u k 2 -2j 2 - Sj u j 2 - 
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Thus the vector S of length 8 1 consisting of the entries 5 tl ^ 2 for the same ordering as 
before is also an eigenvector for A with eigenvalue A = 1. For completeness, we state the 
generalization of Lawton's condition (cf. [10]) in R 2 . 

Theorem 3.1. Let m(x,y) be a given polynomial satisfying ( i) and ( ii) and A a matrix 
defined as in equation ( 59) for the coefficients ofm(x, y). Let (f> be the function generated by 
equation(55). IfX = 1 is a non-degenerate eigenvalue of 'A, then {(f)(-—£i, — ^2) | (^1 , ^2) € 
Z 2 } is an orthonormal set. 



We also need the following well-known necessary and sufficient condition for orthonor- 
mality. 

Theorem 3.2. Let m(x,y) be a given polynomial satisfying (i) and (ii). Let 4> be 
the function generated by equation (55). Then {<p(- — — ti)\[ii,ti) <E Z 2 } is an 
orthonormal set if and only if 

J2 |£((wi,W2) + 27T(fc,*))| 2 = l, V(W1,W2) € [-7T,7T] 2 . 



Case 1: We have used Matlab and Mathematica to check the eigenvalues of the Lawton 
matrix associated with this two-parameter family for a large sample of parameters. The 
eigenvalue A = 1 was non-degenerate for every sample we tested. 

Case 2: These solutions are not associated with scaling functions. The conditions for 
Case 2a immediately imply m(e iwi ,l) = (l+e 5iui )/2.lf we consider the one-dimensional 
restriction fh{uj) := m(e* w , 1), we see that {±-|, ±4f, ±71"} are the zeros of fh(uj). Because 
m(e l ", —1) = 0, condition (ii) implies that |m(cj)| 2 + |m(w + 7r) | 2 = 1. Moreover, 



rn 



3tt 



(-§+-) 



3tt 
I -r 7T 



■rn. 



= 1. 



Because {£1 = ^,£2 = ^,£3 = — X'^ 4 = — T") is a nontrivial cycle in [— tt,u] forthe 
operation ^ — > 2£ mod 2u such that |m(£j)| = 1, the set of functions {4>(- — n)} ne z 
associated with m(w) is not orthonormal and 



£i*( 



2tt 



(See [5]). So, for the unrestricted function, we have 



E 

k,l 



n| y + 2irk, 2tt£ 



= 



which contradicts Theorem 6.2. 

The other subcase 2b has the property that m(e tuJl , 1) = (1 + e t3uJl ) /2 which similarly 
excludes it from being associated with scaling functions. 

Case 3: This case has the same problematic factors as Case 2 but with respect to the 
other component m(l,e laJ2 ) = (l + e i5 " 2 )/2 orm(l,e IU2 ) = (l + e i3 " 2 )/2. 
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Case 4: The solutions for this final case have the same factors as in Cases 2 and 3. 

So, Case 1 is the only solution associated with scaling functions since the other cases 
failed to satisfy the necessary and sufficient condition for othogonality. Although, the 
refinable functions for Cases 2-4 are not orthogonal to their shifts they still have associated 
tight frames since they satisfy condition (ii). These cases are analogous to the univariate 
Haar function with support [0,3]. 



4. THE 8x8 CASE 

In this section, We derive several necessary conditions from the properties (i)-(iv) with 
N = 7 and M = 1. We will use these necessary conditions to show in the next section that 
the second order vanishing moment M = 2 is not possible for this support size. 

7 7 

Let N = 7 and consider m(x,y) = ^ c i jX l y : ' which satisfies properties (i)-(iv) 

i=0 j=0 

with M = 1. The symmetry property (iii) implies that 



m(x,y) 



1 " 


T r 


y 




y 2 




y 3 




y 4 








y 6 
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bo 
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OlO 
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x 4 


b 7 
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bo 
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b 5 
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bi 


a 4 
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_X 7 



Properties (i) and (iii) imply that 



15 1 

E E "< = o- 



«=0 v—a.b 



Property (ii) implies the following 25 nonlinear equations: 

E ^15 = 

i/=a,b 

E ^12 = ° 

v— a,b 

^ I/O ^11 + ^4^15 = 

iy—a,b 

^ ^0^14 + ^15 = 

i/=a,b 

^2^12 + ^3^13 = 

E ^3^8 + ^7^12 = 

E ^0^7 + t^ll + ^8^15 = 



(60) 



(61) 



(62) 
(63) 
(64) 
(65) 
(66) 
(67) 
(68) 



v—a,b 
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^ WW + V\V\i + 2^15 = (69) 
u—a,b 

^ ^1^12 + ^2^13 + ^3^14 = (70) 
u—a,b 

y^ ^3^4 + ^7^8 + ^11^12 = (7i) 

u—a,b 

^ ^0^3 + ^4^7 + ^8^11 + ^12^15 = (72) 

u—a,b 

^ fo^lO + ^1^11 + ^4^14 + ^5^15 = (73) 

^ ^0^12 + ^1^13 + ^2^14 + ^3^15 = (74) 
v— a,b 

^ ^2^8 + 1/3^9 + ^6^12 + ^7^13 = (75) 

i/=a,b 

^ ^0^6 + V\V*J + ^4^10 + ^5^11 + ^8^14 + ^9^15 = (76) 
i/=a,b 

^ V Vg + ^1^10 + ^2^11 + ^13 + ^5^14 + ^6^15 = (77) 
v— a,b 

^ ^1^8 + ^2^9 + I^IO + ^5^12 + ^6^13 + ^7^14 = (78) 
i/=a,b 

^ V2V4 + ^3^5 + ^6^8 + ^7^9 + ^10^12 + ^11^13 = (79) 
v— a,b 

^ VqV S + ViUg + lA^lO + ^3^11 + V4V12 + ^5^13 + ^6^14 
i/=a,b 

= (80) 

y^ vqi>2 + v\V3 + ^4^6 + ^5^7 + ^8^10 + ^9^11 + ^12^14 

v— a, 6 

+^13^15 = (81) 

y^ ^0^5 + ^1^6 + ^2^7 + ^4^9 + ^5^10 + V%V\\ + ^8^13 + ^9^14 
v— a, 6 

+^10/^15 = (82) 

v—a^b 

+viivu = (83) 

^ VQVl + V\V2 + v 2 v z + z/ 4 z/ 5 + is 5 is e + V & V 7 + VgVg + Z/9Z/10 + ^10^11 
v— a, 6 

+2^12^13 + ^13^14 + ^14^15 = (84) 



CONSTRUCTION OF BIVARIATE SYMMETRIC WAVELETS 



29 



^ ^0^4 + ^1^5 + V 2 V§ + V&fij + ^4^8 + ^5^9 + ^6^10 + ^7^11 + ^8^12 
u—a,b 

+VgV 13 + V W V U + ^11^15 = (85) 

EE". 2 = i (86) 

i— u—a^b 

Property (iv) with M = 1 implies 



a + ai + a 2 + a 3 = 6 + h + b 2 + b 3 (87) 

a 4 + a 5 + a 6 + a 7 = b A + b 5 + b 6 + b 7 (88) 

a s + a 9 + oio + an = &s + &9 + &io + &n (89) 

ai2 + ai3 + 014 + ai5 = + h 3 + b 14 + b 15 (90) 

a a + a 4 + a s + a 12 = b 3 + b 7 + b n + b 15 (91) 

ai + a 5 + a g + a i3 = b 2 + b 6 + b w + b 14 (92) 

a 2 + a e + aio + a i4 = b\ + o 5 + 69 + &13 (93) 

«3 + «7 + an + 015 = b + fe 4 + 6 8 + &12. (94) 

Using (87)-(90) and (61), we immediately have 

15 15 



(95) 



i=0 



Next, we use various combinations of the nonlinear equations to make perfect squares 
as we have done previously. This enables us to introduce parameters in order to simplify 
these equations. 

Using (86), (85), (80), and (74), we have 



^ {VQ + V 4 + V$ + V 12 ) 2 + {V\ + V 5 + Vg + Vl 3 y 



v—a.b 



+{v2 + v§ + z/10 + vu) 2 + {u 3 + V7 + v\i + v 15 ) 2 = -. (96) 

Moreover, using (65), (66), (73), (75), (76), (79), (81), we have 
^ (2^0 + v A + v 3 + vi 2 ){v 2 + z/ 6 + v w + v 14 ) 

v—a.b 

+ (^i + ^5 + + vi 3 ){v 3 + vj + vn + fu) = 0. (97) 

After using (91)-(94), (96), and (97), we obtain 

((a + a 4 + a 8 + a 12 ) ± (a 2 + a 6 + a w + a M )) 2 

+((ai + a 5 +a 9 + a 13 ) ± (a 3 + a 7 + an + a 15 )) 2 = -. (98) 
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Choosing the plus sign in equation (98) and using (95) yields 

a + a 4 + a 8 + a 12 + a 2 + a 6 + a w + a\4 = r , (99) 

ai + «5 + «9 + «13 + «3 + «7 + «11 + «15 = SO, (100) 

where r + s = 1/4 and r so = 0. Giving us four cases: r , s — 1/4 or 0. 
Similarly, using (86), (84), (81), and (72), we have 

^ {VQ + Vl+V 2 + V?,) 2 + {"4 + V h + V(, + v 7 f 
u—a.b 

+ (^ 8 + ^9 + ^10 + ^ll) 2 + (^12 + ^13 + ^14 + Vlbf = ^- (101) 

o 

Using (80), (78), (77), (75), (73), (67), (64), we have 

^ (f0 + ^1+^2 + ^3) (^8 + Vo, + V W + l> n ) 
u—a,b 

+{va + V5 + v 6 + v-j){yxi + v lz + u u + v 15 ) = 0. (102) 
In a similar fashion as before, we have 

a + ai + a 2 + a 3 + a$ + a 9 + a w + an = t , (103) 
+ a 5 + a 6 + a 7 + a 12 + a 13 + a u + a 15 = u , (104) 

where to + uq = 1/4 and t^uo = 0. We now refine our solutions to sums of four 
coefficients. Choosing the minus sign in equation (98) yields 

a +a 4 + a 8 + ai2 = -(r + ri), a 2 + a 6 + a w + a u = -(r -ri) 

ai+a 5 + a 9 + ai 3 = -(so + Sij, a 3 + a 7 + a n + a 15 = -(s - Si), 

where ri = | cos a and si = ^ sin a. Furthermore, equations (62), (63), (64), (67), (68), 
(71), and (72) yield 

^(a + a 4 + a 8 + ai 2 )(a 3 + a 7 + an + ai 5 ) = 0. 

This together with (91) and (94), give the following constraint on our parameters 

2(r +ri)(s -«i) = 0. (105) 

Becauseofthe relationship between ri and si, equation (105) produces three cases: ri = r 
and si = s , ri = — r and s\ = s , or n = — r and si = — s - 
Similarly, 

a + ai + a 2 + a 3 = -(to+h), a$ + a 9 + a 10 + an = -(t -ti) 

a 4 + a 5 + a e + a 7 = ^(u + Ui), ai2 + «13 + «14 + ai5 = ^(«o - "l), 

where ti = j cos /3 and ui = -| sin/3. Additionally, equations (62), (63), (65), (69), (70), 
and (74) with (87) and (90), gives us (to + £i)( u o — Wi) = 0. In summary, we have the 
following lemma 
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Lemma 4.1. 

a + a 4 + a 8 + ai2 = 5(7*0 + ri), b + b A + 6 8 + b 12 = 5 ( s - Si), 

ai+a5 + «9 + ai3 = |(so + Si), &i + h + 69 + 613 = \{ro - r\) 

a 2 + a 6 + a 10 + a 14 = \{r - ri), b 2 + b 6 + b w + b 14 = ±(s + si), 

a.3 + a 7 + an + a 15 = ±(s - si), b 3 + b 7 + bn+b 15 = 3 (?"o + »*i), 

a Q + a 1 +a 2 +a 3 = 5 (*0 + ii), b + b ± + b 2 + b 3 = ^(t + ti), 

a A + a 5 + a 6 + a 7 = ±(u + Ui), 64 + 65 + b 6 + b 7 = ±(u + Ui), 

a 8 + a 9 + aio + an = \{t -t{), b$ + b 9 + b w + b n = \{t -t{), 
012 + ai3 + an + ai5 = 5(^0 &12 + &13 + bu + &15 = 5( u o- u i)- 

w/zere 

" , o + s = ^, r s = 0, (r +ri)(s - si) = 0, ri = i cos a, si = ^sina, 
to+uo = ^, t u = 0, (to + *i)(«o - «i) = 0, ti = icos^, ui = j sin/3. 



Because there are no symmetric compactly-supported tensor-product wavelets with more 
than one vanishing moment, it is natural to ask whether there are any nonseparable sym- 
metric solutions with multiple vanishing moments. In [9], it is shown that their are no 
symmetric solutions with higher vanishing moments for m(x, y) with N = 7. Although 
the bivariate case allows enough freedom to generate a family of symmetric solutions, it 
does not allow for multiple vanishing moments at least for the support size we have con- 
sidered. So, compact support, orthogonality, vanishing moments, and symmetry are again 
at odds in the construction of bivariate wavelets. 
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